COOL-DOWN PROBLEM OF A SOLID PROPELLANT 
ROCKET MOTOR INCORPORATING DEWETTING 



Donald Albert Wellmann 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

Cool-down Problem of a Solid Propellant 
Rocket Motor Incorporating Dewetting 

by 

Donald Albert Wellmann 

Thesis Advisor: G. H. Lindsev 



June 1972 



Appfiove.d {^on. pubLLc AeXeoAe; di^tfubixtion uiitarUXzd. 

T147655 



Cool-down Problem of a Solid Propellant 
Rocket Motor Incorporating Dewetting 



by 



Donald Albert Wellmann 
Lieutenant, United States Navy 
BoS., United States Naval Academy, 1964 
M«S», Naval Postgraduate School, 1971 



Submitted in partial fulfillment of the 
requirements for the degree of 



AERONAUTICAL ENGINEER 
from the 

NAVAL POSTGRADUATE SCHOOL 
June 1972 



ABSTRACT 



This thesis undertook a macroscopic, theoretical analysis of a long, 
circular, solid propellant, rocket motor undergoing thermal loading. It 
was shown that classical plasticity's von Mises yield condition adequately 
described dewetting of solid propellant which was modeled to have linear 
strain hardening characteristics. The mechanical properties of the 
propellant permitted a closed-form solution for the stresses and the 
strains resulting from cooling the motor from 160°F to -60^F. 
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I. INTRODUCTION 



The problem analyzed in this thesis was that of a long, cylindrical, 
solid propellant rocket motor (see Figure 1) containing a circular bore 
and with a thin steel case bonded to its outer radius. The motor was 
thermally loaded by cooling it from 160°F to -60°F, which represented 
the extremes that might realistically be expected during its lifetime, 
from cure at the upper temperature to storage at the lower one. The 
cool-down rate was assumed to be slow enough to maintain equilibrium 
conditions throughout the propellant; i. e., temperature did not vary with 
radius. The stresses iii the motor resulted from the difference in the 
coefficients of thermal expansion between the case and the propellant. 
Because the propellant had the higher of the two coefficients, it tended 
to contract more than the case; thus an applied load was induced at the 
outer boundary of the propellant. 

The propellant was considered to be an isotropic, solid material; 
although in reality it was a composite of solid oxidizer particles 
imbedded in and bonded to a carbon fuel matrix. After an initial linear 
region, nonlinear stress-strain behavior under loading resulted from the 
progressive failure of these bonds at the interfaces between the fuel 
matrix and the oxidizer particles. Voids or vacuoles were thus created 
throughout the material, and their volume increased with further loading. 
This phenomenon is called dewetting. 

Dewetting causes nonlinear material behavior and has many similari- 
ties with yielding in classical plasticity theory, as was pointed out 
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by Lindsey and co-workers, (4) who developed an isotropic theory for 
dewettable solids. In it, the strain increment was decomposed into an 
elastic component and a dewetted component. 



d£.: = d£f. + d£;: 



( 1 ) 



•J — ‘i ■ 

They further developed from Drucker’s postulate that the dewetted strain 



increment was normal to the dewetting surface. 









( 2 ) 



where f=o describes the dewetting surface. From this theory. Wood (7) 
showed that the general form of f was 

f= .“.-L (3) 



, I - B for I. ^0 
^ +AI,^ -B for I, > 



where J 2 was the second deviatoric stress invariant, was the first 
stress invariant, and A and B were mechanical properties* Wood experi- 
mentally measured these parameters for a PBAN propellant and found that 
A = -0.0026 and B = 180. 

Lindsey and Wood (A) described the dewetting surface as a circular 

cylinder about the hydrostatic axis in principle stress space. The 

positive end of this cylinder was capped by a body of revolution (see 

2 

Figure 2). Using a space defined by an axis system consisting of 
which was coincident with the hydrostatic axis and which was perpen- 
dicular to the hydrostatic axis. Wood’s dewetting criterion became 
intersecting straight lines with an extremely small negative slope of 
-0.0026 and an intersection with the J 2 axis of 180 , as shown in 
Figure 3. Taking the elastic propellant properties of Wood’s material and 
the elastic solution for an encased cylinder subjected to thermal cool- 
down, the intersection between the thermally induced stress field and 
the dewetting surface was determined by calculating the ratio of J 2 and 
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and finding its intersection with the dewetting surface. 



and I 
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(See Figure 1) It was found that only 1 1/3% error resulted in assuming 
the dewetted surface was a cylinder for the thermal problem. This was 
quite acceptable, especially in view of the great mathematical simpli- 
fication which resulted. Thus, Wood*s dewetting criterion for this 
problem was approximated to be 



where for the thermal problem, B is a function of temperature. It was 
observed that this dewetting criterion was identical to the von Mises 
yield criterion of classical plasticity. 

The goal of this thesis, then was to perform a stress and strain 
analysis using Wood’s dewetting criterion as simplified accounting for 
the nonlinear behavior caused by dewetting, on a long, circular, solid 
propellant rocket motor. 

Prior theory had considered propellant to be completely elastic; 
however, the phenomenon of dewetting was expected .to occur at ‘the 
inner bore initially and to expand outward as loading was increased. 
After yielding occurs in either dewetting or plasticity theory, the 
linear stress-strain law is no longer valid and must be replaced by a 
flow rule. Both theories insist upon a similar form of the flow rule 
and for the dewetting surface of equation (A), the flow law of equation 
(2) becomes ^ 



J2 ^ B(T) 



(A) 




where d£.;^. is the increment of dewetted strain, is the deviatoric 

stress component, andk is a scalar incorporating the loading history. 
Mendelson (5) shows how K can be written in terms of an equivalent 



stress, and o|£pan equivalent dewetting strain increment: 




( 6 ) 
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In this analysis both the dewetting function and the propellant’s 
mechanical properties (except for Poisson ratio which was assumed 
constant) were functions of temperature. These properties were determined 
experimentally for a PBAN propellant manufactured by United Technology 
Center, Sunnyvale, California. Thus in the solution as the temperature 
loading was changed, it was possible to account for variable dewetting 
requirements as well as variable stress-strain relationships. 
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II, BACKGROUND 



Because of the close analogy between dewetting and plasticity, the 
pertinent literature of plasticity was searched for solutions to cylinder 
problems. It was found that the subject of partially plastic, thick- 
walled tubes received much attention from many authors during the late 
forties and the fifties. In general, the vast majority of the work 
dealt with internally pressurized cylinders, which were then of current 
interest and were more amenable to solution than those problems which 
deal with thermal loading. 

In 1957, Steele (6) made a review of several theories for partially 
plastic thick-walled tubes. He stated: 

"The mathematician's approach to the problem under 
consideration is fundamentally different from that of 
the engineer. The former is interested in an accurate 
knowledge of the stresses and strains in terms of the 
applied load or loads for all points in an idealized 
body. In general, the engineer seeks an estimation of 
the variation of factors with applied loads, which 
are of significance only in the functional design of 
the part.... In the case of the thick-walled cylinder, 
interest lies in the deflections of the bore and the 
outer surface for specific internal pressures." 

As an example. Hill, Lee and Tupper (2), in a rigorous work, showed 
that Nadai's approximate solution for the axial stress was in error by 
sixty percent, although the deflections at the bore and outer surface 
correlated closely. Rigorous solutions such as Hill's, Lee's and 
Tupper's, did not yield a closed form solution; and, although mathema- 
tically correct, they were so complicated that their usefulness to the 
engineer was minimal. 
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A general outline of the relationships required by the thick-walled 
cylinder problem is given below. These are required for any mode of 
loading, and apply to elastic, plastic, or dewetting theories equally. 

(1) Radial Equilibrium 

(2) Strain Compatability 

(3) Stress-strain Relation 

(4) Continuity of radial stresses and displacements across the yield 
surface 

(5) Boundary Conditions 

(6) Compressibility 

(7) Yield Criteria 

(8) Yield zone stress-strain relations 

The boundary conditions are fixed by the geometry and the physical 
environment. It was assumed that at the inner bore the radial stress 
was zero and that the outeir surface had an induced pressure caused by 
the restraining effect that the case applied to the propellant as the 
motor v/as cooled. Plane strain was the assumed end condition. 

For mathematical simplicity in many plasticity solutions, the material 
has been assumed to have the same compressibility characteristics in the 
yielded region as in the elastic region. Furthermore, the material 
has been taken to be incompressible, which for metals can introduce 
considerable error. However, this assumption is quite valid when dealing 
with solid propellant, as Wood’s experimental evidence has shown. Thus, 
it was assumed that the material was incompressible in both the dewetted 
and the undewetted regions. This gave considerable mathematical 
simplification in that Poisson’s ratio, 7^, was assumed to be a constant 
and equal to one half and that the sum of the dewetted strain increments 
was equal to zero. 
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( 7 ) 



dt^ dig + d£^ =0 

This assumption permitted the simplification of the yield criterion 
and the stress-strain laws to useable forms. ^ 

Two yield criteria were frequently encountered. The von Mises, 
where 




and the Tresca yield condition, in which 

( 9 ) 

The Tresca condition offered great mathematical simplicity for the 
cylinder problem facilitating solutions that were otherwise out of reach, 
provided that the relative magnitudes of the three principle stresses 
were known a priori. It was shown that in this problem, the axial stress 
was the middle stress, thus there was no need to calculate it. Addi- 
tionally, its contribution to plastic yielding was zero. 

On the other hand, the von Mises criterion required no prior 
knowledge of the stresses, but it yielded nonlinear differential 
equations. Due to the fact that investigators have demonstrated 
experimentally that the von Mises criterion provides a closer correla- 
tion with a description of yield surface than the Tresca condition, 
solutions have been sought and obtained for the pressure problem using 
the von Mises criterion. The thermal problem, however, was more difficult 
in that a thermal term was added to the stress-strain equation. The 
addition prevented a simple closed form solution using the von Mises 
criterion; however, it had no effect on the Tresca condition. Consequently, 
the Tresca condition, even though less accurate, was frequently chosen 
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over the von Mises criterion. Recall, though, that from dewetting 
theory Wood found the dewetting criterion to be of the von Mises form. 
It was therefore desired to obtain a solution using the von Mises 



surface. 

There are two plastic flow rules that have been used exclusively the 
Prandtl Reuss and the Hencky: 

Prandtl-Reuss 



The Prandtl-Reuss equations deal with incremental theory where an 
increment of plastic strain is so small that it can be considered linear. 
The total plastic strain is achieved by integration. An expansion of 
equation (5) for dewetting materials shows it to be identical to the 
Prandtl-Reuss equations . 

Hencky postulated that the total plastic strain would have the same 
type of relationship as the incremental such that the present state of 
plastic strain was independent of its loading path. This is not true 
in general, and the Hencky relationships have fallen into disuse. 

However, under special circumstances of linear strain hardening and a 
monotomically increasing load, the two flow rules coincide (5). 




( 10 ) 



where d£p is the equivalent plastic strain and the equivalent stress. 
Hencky 




( 11 ) 
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Steele presented a summary of the combinations of these parameters that 
have been tried by different researchers. It is repeated here in Table I. 
Note that most of the solutions required numerical analysis. Any 
numerical solution necessitates a loss of generality and makes a para- 
meter study more costly and difficult. Thus a closed form solution was 
desired. By making simplifying approximations which were the result 
of either the special geometry assumed or the unique mechanical properties 
of solid propellant, it was possible to obtain a closed form solution 
to the thermal problem using a von Mises dewetting surface. 
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III. THEORETICAL ANALYSIS 



Within the above mentioned framework and background, the theoretical 
analysis was undertaken. Figure 2 describes the geometry of the problem, 
and the governing equations are given below. 

Equilibrium 

d^r 

dr r ^ 

Strain Compatability 

d £& ^i9~ = 0 

dr r 

Stress-Strain 

e©'- 

where dewetted strain components. 

The problem was well structured in that there were sufficient 
equations to accommodate the unkno\>ms; thus, a unique solution could be 
derived. By applying the plane strain assumption and the incompressi- 
bility condition of the dewetted zone, the last of equations (14) was 
solved for the axial stress 

= l/2( <Dp+<^e)-EOCT 4-E(^ + e‘) (15) 

Using this definition of the axial stress, the first two stress-strain 
relations (equations (14) were substituted into the strain compatability 
equation (equation (13)) yielding a differential equation in terms of 
The tangential stress, was eliminated through the use 
of the equilibrium equation (equation (12)), yielding 



( 12 ) 



(13) 



(14) 
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( 16 ) 



d. i'^'Q - 

cAr i , , 

, 1 dOr _ 2 ( £r-^ )_ „ 
E dr 3 r 



Integrating this equation twice yielded an expression for the radial 
stress. 

^ P P oF I- C .-r’-J.s 1 O 

(17) 



r ^2 3 -4.P 



v‘ ^ IT Jo. 

where and C 2 were constants of integration. Analytic integration is 

only possible when Poisson’s ratio is one-half, which is the actual 

value for propellant. 

The temperature distribution was assumed to be constant as previously 
mentioned. Consequently, the temperature integral of equation (17) was 
readily evaluated. 

^ . 1£l + V 0^ (18) 

r 3 Jo. 'c 

The two constants of integration were evaluated by applying the boundary 
conditions : 

at r = a = 0 

(19) 

at r = b <6r "" 

where p' was the pressure at the boundary between the 

propellant and the case and was caused by the difference between the 

case's coefficient of thermal expansion and that, of the propellant. 

See Appendix A. Application of the first boundary condition yielded 

C. = (20) 

Applying the second boundary condition at r = b, gives 

t-P'’ (21) 

Substituting (19) and (21) back into the expression for yielded 



L ' 3 Jc. ^ 

3 Jc. r 



( 22 ) 
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The need was now to solve the remaining dewetted strain integrals. 



As mentioned above, it had been assumed that the Lindsey dewetted 
stress strain equations applied and that solid propellant possessed 
linear strain hardening characteristics. Expanding equations (5) gives 

d£;.fp(v 

As previously stated, for monotonically increasing loads and a linearly 
strain hardening material, the Prandtl-Reuss equations degenerated to 
the Hencky equations. Within this approximation, therefore, it was 
assumed that the Lindsey equations took the form of 






(24) 



From those equations it is seen that 



where 



t - ^ 



(25) 



(26) 



The axial stress, , was replaced by equation (15), and after expanding 
the terms and regrouping 









(27) 



An examination for order of magnitude of each term of the right side 
of equation (27) demonstrated that the square of the stress difference 
clearly dominated the entire expression. -The radial and tangential 
stresses were of opposite sign and such size that their difference was 
relatively large. The low modulus of the propellant relative to that 
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of some other material such as steel (see Figure 13) made the £=0* term 
small compared to the stress difference. Additionally, from the elastic 
solution, it is seen that the radial strain is negative while the 
tangential strain is positive. Moreover, the absolute magnitude of the 
radial strain is larger than that of the tangential strain; consequently 
the sum of the tangential and radial dewet ted strain components must be 
negative, assuming that the direction of the dewet ted strain increments 
are in the same direction as their elastic components. So it seems 
that the already small &AT term is reduced even more by the addition of 
this negative quantity. Thus, the second term inside the radial was 
neglected. The expression for the equivalent stress became 



( 28 ) 

which simplification made possible a closed form solution without 
appreciable loss in accuracy. The terms neglected were estimated to 
constitute at most a 2% error. 

Using equation (28) , equation (25) reduced to 

(29) 

where for the thermal problem. 

The significance of assuming linear strain hardening actually 
meant that the stress-strain curve was approximated by two straight 
lines, one having the slope of E and the second having the slope of 
mE as shown in Figure 4. From the figure an expression for the dewet ted 
region was found to be 






(30) 



By taking equation (28) and applying the equilibrium equation (12) it 
is seen that 









c\r 



( 31 ) 
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Using equations (31, 30, 29) permitted direct integration of the 
dewetted strain integral in equation (22) such that the radial stress 
expression became 



J _ YY\ 



r_ n'_ ^ 

+ ^5 ^ 4U"/r j 



(32) 



Two new unknown constants appeared in this equation, and r^. 

The critical radius, r^, was the radius of the boundary between the 
dewetted and the elastic zones, and was the radial stress at this 
boundary. Both of these were constant and were functions of the thermal 
loading only. was related to r^ by setting r equal to r^ in equation 
(32). After substituting for p’ 



V 3(\-vtN'^A 



(33) 



where p’ and D came from the elastic solution. (See Appendix A.) 

It was noted that at the boundary between the dewetted and undewetted 
zones, the radial stress must be continuous across the boundary. 

Further, note at that location that the dewetted condition must also be 
met by the stresses in both the elastic and the dewetted regions. 

The yield condition, equation (1) became 

(34) 

The second term under the radical was considered negligible for reasons 
already given. The expression for and were taken from Appendix A. 
The combined loading case was used, and the radius was set equal to the 






critical radius (r=r^) , such that 









(35) 
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Combining equations (34) and (35) and solving for yielded 



- 



r 3 

_ X ^ g.K Uq ^ 



3. _ (l-J4*>bE 
^ e.K 



(36) 



It would have been possible to solve equation (34) and (35) for in 
terms of an unknown r^ and to set this new expression equal to equation 
(33), solving for r^ directly. However, such a course of action would 
have resulted in a transcendental, tenth order expression and would have 
been of little practical value. A simpler method was to iterate with 
equations (36 and (33) . Assume a value for (say = 0) and solve for 
using equation (36). This value for r^ was used in equation (33) 
to solve for a better value of . This process was repeated until 
the resulting values of r^ converged. 

All unknown constants of equation (32) were known and the complete 
expression for the radial stress was obtained: 

(37) 



r>rt 






(38) 



Through the application of the equilibrium equation, the tangential 
stress was obtained directly. 






(39) 



The strains were obtained by application of the stress-strain 
equations. 



c _ _ 3 _ (<^r ~Q ^ ^ d(T _ ^ 

M VY\ E Z. HExt\ 

- •- ^ 



(40) 



e 



Vt\ 



MErrx 
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IV. EXPERIMENTAL RESULTS 



There was a dual goal in the experimental work done. First, it 
was necessary to experimentally determine the constant, B(T) in equation 
(2), for the particular PBAN propellant used. This was accomplished by 
conducting uniaxial tensile tests over a range of temperatures in the 
gas dilatometer of United Technology Center, S.unnyvale, California. 



For uniaxial tension. Equation (2) reduced to: 
B(T) = % 



(41) 



The dilatometer (see Figures 5, 6, 7, 8) utilized the Perfect Gas 
Law as its basic operational principle. The specimen was mounted within 
a sealed chamber and pull-rods extended from either end of the sample 
through the ends of the chamber and attached to the cross-head of an 
Instron testing machine. Any change in the volume of the sample while 
under load resulted in a change in the chamber’s pressure creating a 
pressure difference across a diaphram separating the test chamber from 
an adjoining reference cavity. A baratron pressure transducer converted 
this pressure difference into an electrical signal, which was fed to 
one pen of a two-pen X-Y plotter. The signal for the second pen came 
from a load cell in the cross-head of the Instron machine. The raw data 
thus obtained were in the form of two superimposed plots; one was the time 
history of the load, and the second was the time history of the dilata- 
tion. By using the constant cross-head speed and premeasuring both 
the cross sectional area and gage length, the data was reduced to plots 
of stress vs strain and dilatation vs strain. (See Figure 9) 
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Changes in the ambient temperature caused a problem since it 
affected the readings of the pressure transducer through the same Perfect 
Gas Law. In order to isolate the experiment from ambient temperature 
fluctuations and also provide a method of controlling the temperature 
at which the experiment was conducted, the dilatometer was enclosed in 
a temperature box, which possessed the capability of maintaining a 
preselected temperature to within a degree over the duration of the 
experiment. Although the ambient temperature within the temperature 
box was fairly constant, at subzero temperatures, heat leaked into the 
dilatometer probably by conduction through the push rods. A large 
temperature drift would develop after about one half hour and this 
dictated the maximum duration of any individual run. 

The second goal was to experimentally determine the mechanical 
properties of the subject propellant as functions of temperature; 
namely, the uniaxial dewetting stress, , and the elastic modules, E 
and m. (See Figures 10, 11, 12) 

The test samples were square with sides machined to approximately 
one half inch, and they had a gage length of approximately 23/4 inches. 
Redwood tabs were bonded to each end of the sample to make hardware 
attachment possible. Tests were conducted at 21°C, 3®C, -9^C, -21^C, 
and -52®C, with four tests being conducted at the three higher tempera- 
tures, two tests conducted at -21® C and one test conducted at -52°C. 
Experimental difficulties with the equipment at the subzero temperatures 
resulted in the reduced number of runs at the lower temperature. 
Experimental scatter of seven to ten percent was found in the dilatation 
and stress-strain curves. However, only about five percent scatter was 
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found among the parameters measured (i.e., E, m, and 6^ ) . Due to the 
fragile nature of solid propellant specimens, this was considered quite 
reasonable results. 

After reducing the raw data to the stress vs strain and dilatation 
vs strain plots mentioned above, the problem of identifying developed . 
Three methods were available from the literature: (1) interpreting the 

stress which corresponds to the lift-off point of the dilatation curve 
as . This method provided results which were quite repeatable; 

however, the stress strain curve did not depart from its initial linear 
relationship at this point. (2) The dilatation curves eventually 
assume a constant slope. Physically the total amount of solid material 
which could dewet had dewetted at this point and further dilatation 
resulted only from increasing the size of previously formed voids. The 
stress corresponding to the extrapolation of this constant slope line 
back to the zero stress point (pt. B), of Figure 13 would be considred as 
the yield stress, ^ . This could be physically interpreted as assuming 

o 

that no dewetting occurred until the yield stress was reached; then all 
the dewetting occurred simultaneously. This was tantamount to assuming 
that the propellant was perfectly plastic. However, at the lower 
temperatures, specimen failure resulted before such a linear portion to 
the dilatation curve was reached. (3) The third method v;as to take the 
stress at the point where the stress-strain curve initially departed 
from a linear relationship. This method reflected the linear-strain 
hardening assumption, and provided consistent reproduceable results 
from sample to sample. It neglected the dilatation curve altogether. 
This last method was adopted primarily because it met the condition 
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of the linear strain hardening assumption. A comparison of these 
three methods was given in Table II, where it is seen that methods 
(2) and (3) yield similar results. 

The test speeds were chosen as slow as possible, yet fracture had to 
be reached within the time limit mentioned above. At the lower tempera- 
tures, it was assumed that the speeds used maintained the specimens in 
equilibrium and that the load did not produce viscoelastic effects. 

In summary, the yield stress, <5^ , was taken as the stress at which 
the average stress strain curve first departed from linearity for all 
runs at any one temperature. The modulus, E, was the average of the 
elastic slopes of the stress-strain diagrams for any given temperature. 
The strain hardening exponent, m, was the ratio of the slope for the 
least squares straight line fit of the non-linear portion of the stress 
strain curve, to the initial modulus. 
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V. RESULTS AND CONCLUSIONS 



A. RESULTS 

In order to compare the results predicted by dewetting theory with 
those obtained from linear elastic theory, a sample problem was worked. 

A PBAN propellant of current interest was used, and its mechanical 
properties were experimentally determined, (Figures 11, 12). The other 
required parameters were taken from Appendix II of Reference (4) . 

For the propellant: ~U - 1/2 ck- 2.21x10"“^ in./in./®F 

For the steel case: = 1/3 .65x10"*^ in./inc./®F 

E = 3xl0^psi 
c 

Geometrical parameters: a= 3in. b= 12in. h=0.015in. 

The rocket motor was cooled from 160°F to -*60°F, and both elastic 
and dewetted stresses train analyses were made. Linear strain hardening 
was assumed which enabled calculations to be performed in two steps: 
one elastic step and one dewetted step. The temperature was not 
incrementally stepped down. Equations (32) and (36) were iterated in 
order to find the critical radius, r , and the interface stress, <^r . 
These results were plotted in Figure (14). It was interesting to note 
that the growth of the yield surface was linear with temperature. 
Furthermore, dewetting did not begin until the temperature had been 
decreased to 10® F, and at the lower temperature limit the yield surface 
had penetrated to a radius of only 9 1/2 inches. The interface pressure 
also increased linearly as the temperature was lowered reflecting the 
influence of temperature on both the yield criterion and the propellant’s 
mechanical properties. 
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The greatest effect of dewetting was upon the tangential and radial 
stresses. (Figures 15, 16) The tangential stress at the bore, for both 
the elastic and the dewetting theories, was plotted against temperature. 
It was seen that dewetting substantially relieved the stress. The 
radial stress was plotted across the motor ^s thickness for a characteris- 
tic temperature of -8.0®F. Again it was shown that stress was relieved 
within the dewet ted region. The radial stress then approached the 
elastic solution as one moved deeper into the undewetted elastic 
material and away from the critical radius. However, a pronounced 
effect of the induced pressure at the inner radius of the undewetted 
region was still evident. 

Initially, the behavior of the tangential strain at the bore was 
quite surprising. (See Figure (17)) It had been anticipated that this 
strain would have been negative and would have increased in magnitude 
as the temperature was lowered. How^ever upon looking at the elastic 
tangential strain equation (A-1), one sees that it consists of two terms 
The first term represents the mechanical strain imposed by the case 
upon the propellant. It is positive and dominates the second term, 
the thermal strain, which is negative in the cool-down problem. However 
as the temperature is lowered the mechanical properties of the pro- 
pellant change significantly, and the negative T term grows in relative 
magnitude, actually forcing the total strain to decrease. Thus, it was 
observed that the maximum strain occurred at about 40®F instead of at 
the lower temperature extremum. 
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B, CONCLUSIONS 



A point worthy of note is the curve fit of the m curve coefficient 
of linear strain hardening in figure 12. Equipment limitations at the 
subzero temperatures, which restricted multiple runs decreased confidence 
in the values of m. Consequently, the specific mechanical properties 
measured at -5.8°F and -60°F were used to calculate the bore hoop 
strain at those two temperatures, and they appear on the plot in Figure 14 
which uses the curve fit for m. It was concluded that a curve fit which 
averaged these values provided, a smoother extension of the elastic 
solution than a curve which connected each data point. Further tests 
would prove useful to more definitely fix the m curve. 

From a comparison of the elastic and the dewetting solutions for 
both stresses and strains it was concluded that dewetting has little 
effect on the strains. In an incompressible material the strains are 
fixed by geometrical constraints and not by mechanical properties. 

However, stresses are tied to the mechanical properties and thus are 

greatly affected by dewetting. 

« 

It was further concluded that a thermo-dewetting solution could 
be obtained in a relatively simple form and that better results for stress 
would be obtained by using dewetting theory vice elastic theory. 
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COMPARISON OF SOLUTIONS FOR THE THICK-WALLED PLASTIC CYLINDER PROBLEM 
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° Inclusion of strain-hardening; nonas terisked solutions assume material to flow according to function 
defining theory of failure. 

^Solution dependent on experimental measurements which Cook made on closed-ended cylinders. 
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TABLE II 



TABULATION OF EXPERIMENTAL DATA 



Temp 


cr 

o 


E 


m 


21° 


94 


806 


,417 


3.5 


118 


983 


.545 


-9 


137 


1600 


,49 


-21 


159 


8840 


,369 


-52,5 


231 


23100 


.539 
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TABLE III. 



COMPARISON OF METHODS OF READING 
TENSILE YIELD STRESS 



TEMPERATURE 

(°c) 


I 


^ (psi) 
II 


III 


22° 


46 


84.5 


94 


3.5 


72.5 


107 


118 


-9 


105 


★ 


137 


-21 


137 


★ 


159 


-52.5 


201 


★ 


231 



^Dilatation Curve did not reach constant slope. 



I® <^o was stress corresponding to the lift off point 
of the Dilatation Curve. 

II® 6c was stress corresponding to extrapolation of 

of the constant Dilatation Curve back to the strain 
axis ® 

III® 6c was stress at which stress-strain curve initially 
departed from a linear relationship. 
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Dewetted Propellant 



Rocket Motor Geometry 
Figure 1 . 
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Dewet ted Yield Surface in Principle Stress Space 

Figure 2, . 
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Intersection of Thermal Stress Load 
with Dewetting Surface 

Figure 3, 
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Linear Strain Hardening 
Figure 4. 
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Figure 5 
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Figure 



Figure 7 
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Figure 
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UNIAXIAL TENSILE STRESS VS STRAIN A>ILATATI0N VS STRAIN 

(Temp = 20°F) 
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UNIAXIAL TENSILE STRESS VS STRAIN 
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UNIAXIAL TENSILE YIELD STRESS VS. TEMPERATURE 
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MODULUS AND LINEAR STRAIN HARDENING COEFFICIENT VS TEMPERATURE 
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Calculated Uniaxial Tensile Yield Stress 
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Critical Radius and Interface Pressure vs Temperature 



60(V 
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Figure 15 . 
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Radial Stress 



Dewetted Solution 
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Tangential Strain vs Temperature 
Figure 17. 



APPENDIX A 



ELASTIC SOLUTION 



Williams, Blatz and Schapery (7) of the Guggenheim Aeronautical 
Laboratory (California Institute of Technology) analyzed solid pro- 
pellant rocket motors of the same geometry under many different loading 
modes. Their basic assumptions were that 

1. the propellant had material properties x^hich were 

a. Isotropic 

b . Homogeneous 
d. Continuous 

2. the strains were small. 

Their solution to the thermally loaded, cylindrical, rocket motor, 
assuming a condition of plane strain, was: 




(A-1) 



where 
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I 



And for the case of an internally pressurized rocket motor their 
solution was ; 






^ - ficOL-Pcb^ 

0 \t y^-L 



3 q} 4>i 

2 “CL> 



(A-2) 



D 



where = internal pressure 



In considering the effects of dewetting, the inner boundary condition 
changed. The inner radius of the elastic region was the same as the 
critical radius, i.e., it was equal to the radius of the dewetting 
surface. As the dewetting surface* moved outward, the inner radius 
of the elastic zone also increased; furthermore the inner boundary was 
no longer stress free. From continuity of the radial stress across the 
dewetting surface, it was seen that the inner boundary now had a stress 
equal to ^r^ induced upon it. This elastic solution was obtained by com- 
bining A-1 and A-2, and setting p^ equal toC^r^, and the inner radius, 
a, equal to the critical radius, r^ 
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